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Abstract 



We calculate the exact Green function of a special model of coupled Lut- 
tinger chains with arbitrary interchain hopping t±. The model is exactly 
solvable via bosonization if the interchain interaction does not fall off in the 
direction perpendicular to the chains. For any finite N we find Luttinger liq- 
uid behavior and explicitly calculate the anomalous dimension 7^^^. However, 
the Luttinger liquid state does not preclude coherent interchain hopping. We 
also show that ^^^^ — > for ^ 00, so that in the limit of infinitely many 
chains we obtain a Fermi liquid. 

PACS numbers: 71.27.-ha, 05.30Fk, 71.20.-b, 79.60.-i 
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Typeset using REVTeX 



The normal metallic state of interacting electrons in one spatial dimension is called a 
Luttinger liquid, and has fundamentally different properties than the Fermi liquid state in 
higher-dimensional Fermi systems. In particular, the single-particle Green function of a Lut- 
tinger liquid does not have a quasi-particle peak, but exhibits spin-charge separation and 
anomalous scaling properties. In order to calculate the physical properties of Luttinger liq- 
uids, non-perturbative methods are necessary, such as the bosonization method . Recently 
the crossover from Luttinger liquid to Fermi liquid behavior as function of some suitable 
parameter has received much attention p|, partially because of Anderson's and coworkers' 
suggestion that higher- dimensional Luttinger liquid behavior might be the key to under- 
stand high-temperature superconductivity. An experimentally important class of systems 
exhibiting the above mentioned crossover are quasi-one- dimensional chain-like metals. In 
the normal metallic state (which is realized at sufficiently high temperatures), these systems 
can be described as one-dimensional Luttinger liquids that are coupled by small interchain 
hopping t± and by the three dimensional interaction. While for = the number of parti- 
cles on each chain is conserved and the problem can be studied in a straightforward way via 
bosonization, the case of t_L 7^ is more difficult and has recently been studied by means 
of a variety of non-perturbative methods. Most authors have focused on the crossover as a 
function of the dimensionless parameter t±/Ep (where Ep is the Fermi energy), assuming 
that this parameter is small. In this work we shall show that for some special type of in- 
teractions (see below) it is possible to obtain an exact solution for the Green function of an 
arbitrary number of coupled Luttinger chains with finite t± . The interaction is characterized 
by the fact that it does not fall off in the transverse direction. Although such an interaction 
is unphysical, our model is useful for testing the various approximations employed in the 
literature. 

Let us start with a two-dimensional system of interacting fermions moving on an array 
of weakly coupled spin 1/2 chains in the x-direction. Each chain can be described by the 
Tomonaga-Luttinger model with spin [Q. Denoting by k\\ the momentum along the chains 
and linearizing the energy dispersion in the chain direction, the kinetic energy part of the 
Hamiltonian can be written as 

^kin = EEE [kj^'^ih) - M {^^ih)^Uh) - {£?Kh)^Uh))o) , (i) 

where c°/(/c||) and cf^{k^\) are the creation and annihilation operators for the right and left 
moving fermions (a = +, — ) on chain i = [— y + l], . . . , [y] with spin a =1, | and momentum 
k\\ in the chain direction. The linearized energy dispersions are e^(/c||) = VF{±k — kp), where 
vp is the Fermi velocity and kp is the Fermi momentum associated in the absence of inter- 
chain hopping. For simplicity we assume only nearest neighbor hopping and impose periodic 
boundary conditions in the transverse direction. In this case the interchain hopping is for 
A^ > 2 of the form tij = tj-j,o = <^|i-i|,i^±/2. In the special case of just two chains periodic 
boundary conditions are automatically satisfied and we set tij = t±S\i-j\^i. Because of the dis- 
crete translational invariance in the transverse direction (the y-direction) we can completely 
diagonalize i^kin via a discrete Fourier transformation in the transverse direction. Denoting 
by k± the corresponding transverse momentum and defining c"(k) = A^~^/^ c~*'^^^'c"g.(A;||), 
we have 
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(2) 



where k = (fcy, k±), and the k-dependent energy dispersion is given by e"(k) = e°'{k\\)—t{k±), 
with t{kj_) = tj_ cos{k±aj_). Here a_L is the distance between the chains. The corresponding 
Fermi surface is shown in Fig.^ 

The total Hamihonian is given by if = if km + -f^mt? where the interaction part is in 
general of the form 



1 1 



(3) 



i,j a,a' a,cr' g\\ 



Here, one factor of 1/2 has been introduced for convenience due to the spin, and the Fourier 
components of the density operators associated with chain i are given by 



ZW) = E i^^ikKAk + q)- {^HkKAk + q) 



(4) 



Assuming spin and inversion symmetry, we thus assume that the interaction parameters 



satisfy /+^t = f„„, , = /^J, f{\ 



TT 



the complete Fourier transformation of iiint by defining 

= ^ E (c"nk)c^(k + q) - (c^t(k)a«(k 



/ff', and /^"' = /ff' = /ff'. Performing 

q))o) , 



k 

/r'(q) = Ee-^'"'"/XU?ii 

i 

and going over to a charge/spin basis. 



(5) 
(6) 



Pc/s(q) = ^K(q)±pr(q) 

/c"/:'(q) = J[/ir'(q)±/r'(q) 



we obtain 



iiir 
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E E E/r'(q):pr(qK'(q): 



(7) 
(8) 

(9) 



a,a' v=c,s q 



We would like to calculate the single-particle Green function of the model defined above. 
In particular, we are interested in the fate of the Luttinger liquid state as function of t± 
and N . Note that for t_L = the above model is exactly solvable by means of bosonization, 
because in this case the particle number on each chain is conserved [Q. The case of finite 
is very difficult to handle, and there exists no complete agreement in the literature about 
the nature of the ground state. In this work we would like to point out that there exists a 
special type of interaction (q) where the above model is exactly solvable via bosoniza- 
tion for arbitrary tj_ and N. This is easily seen in the functional bosonization approach 
0J^, where the interaction part of the effective action corresponding to ifint is decoupled 
via a dynamic Hubbard-Stratonovich field (p. In this approach the exact solubility of the 
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Tomonaga-Luttinger model manifests itself via the fact that in the perturbative expansion 
of the effective action for the 0-field, which is obtained by integration over the fermionic 
degrees of freedom in the usual way, all non-Gaussian terms vanish identically (see Fig. 
H). This cancellation has first been noticed by Dzyaloshinskii and Larkin 0, and was later 
discussed in detail by T. Bohr p], who formulated this calculation in terms of a theorem 
which he called closed loop theorem. While this theorem is exact in the one-dimensional 
Tomonaga-Luttinger model, it remains approximately valid even in higher dimensions if the 
interaction is dominated by forward scattering 

The crucial observation of this work is that for our model defined above the closed 
loop theorem is still exact for interactions of the type (q) = for qx ^ 0, which in 
real space amounts to a potential which is independent of the chain indices. In this case 
the auxiliary fields do not transfer any transverse momentum into closed fermion loops. 
For a linearized energy dispersion along the chain direction, the loops with more than two 
external fields cancel then for exactly the same reason as in the one- dimensional Tomonaga- 
Luttinger model. After Fourier transformation this system is essentially equivalent to a 
system consisting of N independent Luttinger chains each with a different Fermi-momentum 
and with only one Luttinger chain showing interaction. We would like to emphasize that this 
model is exactly solvable for arbitrary and t^, although it is physically meaningful only 
for \t_[_\ -C Ep, because we have linearized the energy dispersion along the chain direction. 

The calculation of the imaginary-time Green function is now analogous to the calculation 
of the Green function of the Tomonaga-Luttinger model. For simplicity we set the chemical 
potential /i = 0. In this case there is no need for a special treatment of the particle mode 
p!0| . Following a suggestion of Luther and Peschel |1T[ , we assume a potential which satisfies 
= 7i/e~^'''ii', where 



vZ\^) =vf + 



^^(g|hO) _ ^ 



27r 



«(q))^ - 



• J V 



2-K 



(10) 

(11) 
(12) 



For L 
tures 



oo we obtain for the real-space imaginary-time Green function at finite tempera- 



*(fcx). 



exp (Q;^(x,r)) 



(13) 



where 



a 



{7r/pVF)e 



iakpx 



exp (Q°(a;,r)) = L„{x,t)K„{x,t 



2'Ki smh[{T: / j3vF){x + iavpT)] 
X + iavpT + iar 



X + iaVyT + iar 



{x + iVyT + ir){x — iVyT — ir) 



(14) 



(15) 
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LJx.t) 



T{aF + iup)T{aF — iup) 
r(a^ + w+)r(at, - iu^)r{a^ + iu~)r{a^ - iu~) 
V{ay)V{ay)V{ay)V(a^) 



(16) 
(17) 



and ttj = 1 + r(3vi^uf = (x + iavir) / {(3vi). Note that at zero temperature both factors 

L^{x,t) and K^{x,t) are equal to one. Because the only fc^-dependence is in the e*" ^- 
factor, this Green function is easily transformed to the real-space basis in the transverse 
direction, 



gf,iJx,r) = gS{x,T)g,^i{x) n 



exp (Q;^(x,r)) 



1 

2N 



with 



9j-iix) 



1 \ ik±a±{j ~l)+ia^^x cos{k±aj_) 
k I 



(18) 



(19) 



This is the exact expression for the Green function in real space and imaginary time. The 
only )f:_|_-dependence is in the factor gj_i{x). In the special case of just two chains this factor 
simplifies for j = I to cos and for j ^ I to (ia) sin (^^x^. 

The above result has a number of interesting properties. For any finite the Green 
function in Eq. ([T8|) shows typical Luttinger liquid behavior. There are different possibilities 
of scaling the interaction as the number of chains increases. Analogous to the Weiss-model 

one could choose f°j" (gy) = " where (gy) is independent of A^. In this case 



/""'(q) = ^q±,ofu"' (Qw) is independent of A^, and for A^ ^ oo the anomalous dimension 
trivially scales to zero. Then the Green function of the interacting model reduces to the 
Green function of the noninteracting model. We will now show that this is even true for 
a model for which the interaction between two given chains does not depend on the total 
number of chains A^, i.e. (gy) = (gy). In this case we have (q) = Sq^^N f""" (gy). 
The anomalous dimensions of the charge and spin channel can then be read off from Eq. 
(IMID and Eq. (|18D, 



7^^ 



2N 



vf + N 



/.++(o) 



2n 



Vf 



(20) 



For N ^ oo the anomalous dimensions still vanish. For = we have 
while 7^^^ oc A^~^ for /++ > fl^~. In fact, it is easy to show that in both cases 



lim [exp{Q'^{x,T))]N = 1 

N^co 



(21) 



so that in this limit the electron-electron interaction does not affect the Green function at 
all. Then Eq. (p!8| ) reduces to the exact non-interacting result with hopping, 
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lim 

N^oc 



gU^ix,T) =^-(a;,r)Ml^-'lj|,_,|(^x) , (22) 

Vp 



where Jp (a) is the Bessel function of first kind and order p. Obviously, we are not expanding 
in powers of t±. Thus, while for any finite N the Green function shows Luttinger liquid 
behavior, for iV — >^ cxo it reduces to the Green function of a trivial Fermi liquid. 

To understand this at the first sight rather surprising result, consider the leading self- 
energy correction E"(k, in an expansion in powers of the RPA (random-phase approx- 
imation) interaction. Assuming for simplicity that is independent of a and a' and 
omitting these indices, we have 

S-(k, = Y: ;7^4^^o (k - q, - lu^m) , (23) 

P^^-^ (\,ujm ^^rn) 

where Un is a fermionic and Um a bosonic Matsubara frequency, Gq (k? i^n) = [i<^n — ^"(X)]~^, 
and the RPA dielectric function is given by 

e(q,z..J = l + ^^4^ . (24) 

Note that the closed loop theorem guarantees that all corrections to the dielectric function 
beyond the RPA cancel. Our model corresponds to an interaction /^(q) = 6q^flNf^{q\\). 
Hence only the gy-summation in Eq. (p3D survives, and the prefactor of is canceled by 
the factor N in the bare interaction. The above self-energy reduces then to the self-energy 
of the Tomonaga-Luttinger model, except that the effective dielectric function diverges in the 
limit N oo. Thus, while for any finite N the perturbative self-energy of our model exhibits 
the same type of singularities as the Tomonaga-Luttinger model, in the limit N —>■ (yo the 
effective interaction is completely screened, so that the lowest order self-energy vanishes. 
Our exact solution shows that this is true to all orders in perturbation theory. 

The strong screening of the interaction is clearly a consequence of the fact that our model 
has a super-long-range bare interaction that does not fall off in the transverse direction. Al- 
though interactions of this type are unphysical, we have learned two conceptually important 
points from our calculation. First of all, a system consisting of a finite number of chains can 
exhibit qualitatively different behavior from a system of infinitely many chains. In partic- 
ular, conclusions drawn from two coupled Luttinger chains are in general not applicable to 
the physical more relevant case of infinitely many chains. In our model, we obtain Luttinger 
liquid behavior for any finite A^, but a trivial Fermi liquid for — > oo. As shown above, 
for large A^, the anomalous dimension vanishes as N~^^'^ for = /^~, or as A^~^ for 
> In numerical investigations this scaling of the anomalous dimension with the 

number A^ of coupled chains could be a useful guide for the extrapolation to the infinite chain 
limit. Finally, we would like to point out that two exactly solvable two-chain models that 
have independently been proposed by Shannon, Li and d'Ambrumenil [|T^| can be obtained 



as special cases of our more general model with two chains (A^ = 2) by choosing special 
types of interactions. As first noticed by Shannon et ai, these models have the remarkable 
property of exhibiting Luttinger liquid behavior together with coherent intrachain hopping. 
This is also true for our more general model of A^ coupled Luttinger chains with spin, as 
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is evident from the factorized form of the Green function given in Eq. (|T8|). Thus, it is in 
general not correct that Luttinger hquid behavior in an array of coupled chains precludes 
coherent single-particle hopping . However, this might be a special feature of our model, 
which has the for a large number of chains unrealistic property that the interaction does 
not fall off in the transverse direction. As shown here, the two-chain models proposed by 
Shannon et al belong to the general class of models of coupled chains that are character- 
ized by an interaction which does not transfer any transverse momentum. These models 
are equivalent to a system of N independent one-dimensional chains. We believe that the 
coexistence of coherent interchain hopping and Luttinger-liquid behavior is a special feature 
of these models and is not relevant to physically more realistic models where the interaction 
has a finite range in the transverse direction . 

We would like to thank Kurt Schonhammer for helpful discussions. 
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FIGURES 



FIG. 1. Fermi-surface for the model with N = 6 chains and periodic boundary conditions. 




FIG. 2. Feynman diagram for a closed loop. The solid lines with an arrow represent free Green 
functions and the wiggled lines represent the auxiliary field 4>. The closed loop theorem states that 
the contributions from all diagrams with more than two auxiliary fields in a closed loop vanish. 
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